Slonczewski windmill with dissipation and asymmetry 
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J. Slonczewski invented spin-transfer effect in layered systems in 1996. Among his first predictions 
was the regime of "windmill motion" of a perfectly symmetric spin valve where the magnetizations 
of the layers rotate in a fixed plane keeping the angle between them constant. Since "windmill" 
was predicted to happen in the case of zero magnetic anisotropy, while in most experimental setups 
the anisotropy is significant, the phenomenon was not a subject of much research. However, the 
behavior of the magnetically isotropic device is related to the interesting question of current induced 
ferromagnetism and is worth more attention. Here we study the windmill regime in the presence 
of dissipation, exchange interaction, and layer asymmetry. It is shown that the windmill rotation 
is almost always destroyed by those effects, except for a single special value of electric current, 
determined by the parameters of the device. 

PACS numbers: 72.25.-b, 85.75.-d 



Spin-transfer effect as a method of controlling mag- 
netic dynamics by electric current was suggested by 
Berger— for domain wall motion and by Slonczewski^ for 
spin-valves and multilayer structures. The unusual prop- 
erty of spin-transfer interaction found in Ref. was the 
tendency of current induced torques to rotate magnetic 
moments of both spin valve layers in the same direction, 
much like the oncoming wind rotates the wings of a wind- 
mill (Fig.[T|). If one assumes that layers have no magnetic 
anisotropy (crystalline or shape), are identical, and there 
is no RKKY exchange or dipole-dipole interaction be- 
tween them, the resulting motion is a perpetual rotation 
of magnetic moments mi and m2 in clockwise or counter- 
clockwise direction, depending on the direction of electric 
current I passing through the spin valve. We will call this 
type of motion a Slonczewski "windmill regime" . 

Since actual spin-transfer devices have significant mag- 
netic shape anisotropy, normally the windmill regime is 
not realized. Instead, switching between different pre- 
ferred magnetic configurations was predicted^ and is in- 
tensively studied since then both experimentally and the- 
oretically. However, the windmill regime still constitutes 
an interesting problem due to the following. Spin-transfer 
effect can be viewed as reciprocal to the giant magne- 




FIG. 1: (A) spin-transfer torques acting on two magnetiza- 
tions tend to rotate them in the same direction. Magnetic 
layers are shaded. (B) "windmill" mechanical analogy. 



toresistance effect."^ The resistance of a spin valve is min- 
imal in the parallel configuration mi ]] m2. Thus one 
can hypothesize, that in response to a current pumped 
through the valve the magnetizations will tend to as- 
sume this minimal resistance configuration in order to 
make electron flow easier. More generally, an idea arises 
that a current passing through a metal with paramag- 
netic impurities will tend to orient them parallel and cre- 
ate some sort of current-induced ferromagnetism^i? The 
two-magnet device is the minimal model where the valid- 
ity of this idea can be tested. We study the behavior of 
such a device with arbitrary parameters, except for the 
restriction of zero magnetic anisotropy. The results give 
a generalized picture of the Slonczewski windmill regime, 
and shed some light on the possibility of current-induced 
ferromagnetism. 

We use the single domain approximation. The mag- 
netic moments m^ (i = 1,2) of the layers have time- 
independent absolute values rrii and variable directions 
defined by a unit vector nj(i). The LLG equations in 
terms of m.^ read^ii^ 

Oil 

rhi = 7(Tea;+-ri)H [mi x rhi] , (1) 

mi 

m2 = l{-Tex + T2)^ [m2 X 012] , (2) 

m2 

where Tex is the exchange torque, ti^2 are spin-transfer 
torques, 7 is the gyromagnetic ratio, and ai_2 are Gilbert 
damping constants of the magnets. Note that in conven- 
tional experiments m2 is fixed by magnetic anisotropy, 
while mi can rotate under the influence of spin-transfer 
torque. Magnet number one is then called a "free layer" 
and magnet number two is called a "fixed layer" , or spin 
polarizer. In the present investigation no restrictions are 
imposed on m2 and both magnetic moments are treated 
on equal footing. 

The exchange torque acting on mi is given by Tex = 
J[m2 xmi] ( J > corresponds to ferromagnetic coupling 
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between the moments). The exchange torque acting on 
m2 is —Tex since we are deaUng with an internal inter- 
action between two moments. 

The spin-transfer torques T12 are given by 



Ti = Ml • [ni X [n2 X niJJ , 
T2 = -U2 ■ [n2 X [ni X n2] 



with torque strengths 
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g,[(ni • n2)] 



(3) 
(4) 
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Here / is the electric current flowing form magnet 
2 to magnet 1, e is the (negative) electron charge, 
and gilini ■ 112)] are material and device specific spin- 
polarization factors. For negligible spin-relaxation in 
the non-magnetic spacer between the magnets one has 
gi = 92 (see Ref. 121). Note that both ui,2 are positive 
when electrons flow from magnet 2 to magnet 1. The mi- 
nus sign in front of the right hand side of Eq. Q reflects 
the symmetry of spin-transfer torque.^ 

First, we rewrite Eqs. ^ and ^ so that time deriva- 
tives are on the left hand side only. Deflning Ti = 



Ti, 



T2, we get for i = 1,2 



(1 + af)m, = 7 T, + — [mi X T,] 



(6) 



It is convenient to introduce vectors = [ni2 x mi], li 
[mi x [m2 X mi]], I2 = [m2 x [mi x m2]] . Then 

rill = Aiv + Bill , 
rii2 = -A2U + B2I2 , 



(7) 



with 



Ai = 
A2 = 
Bi = 
B2 = 



l + al 

7 

l + al 

1 



7 I J aiui 



J 



(1 + a\)mi 



mim2 
0:2^2 
mim2 

aiJ + 



Ui 



r71lTO2 
U2 



\ Ci2J 

(1 + a^)m2 \ mim2 
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Consider now the total magnetic moment of the system 
M = mi + m2 and calculate the derivative dWP /dt. Us- 
ing Eq. ([6]) and the properties (m^ • 1^) = 0, ((mi -l-m2) • 
li) =v^, wc flnd 



dt 



= 2C[mi X ma]^ 



(9) 



with constant coefficient C that depends on material pa- 
rameters and spin-transfer strengths 



C = Bi + B2=i 



ai 



a2 



(1 + q;5;)toi (l + a2)TO2 

1 I Ui U2 



mim2 \ (1 -I- ai)TOi (I + al)m2 



J 



(10) 



Since [mi x m2]^ is always positive, except in parallel 
or antiparallel configurations, we can conclude that after 
a transient period the magnetic configuration will reach 
either the state of maximal M (i.e., parallel state) for 
C > 0, or the state of minimal M (i.e, antiparallel state) 
for C < 0. Since in both coUinear states Tex — and 
T1.2 — 0, the system will come to rest and no "windmill" 
motion will happen. For small spin transfer torques ui^2 
the final state will be determined by the sign of J and, 
as expected, the device will end up in a configuration 
corresponding to the minimum of exchange energy. 

The marginal case C = is the only situation when 
the "windmill" is possible. According to Eq. ([9]), the 
value of C linearly depends on electric current / through 
Ml. 2- The only exception is the singular case when device 
parameters satisfy 5i/[(l + Q;i)mi] = g2/[(l + a|)m2], 
and C is current-independent. Thus in general one can 
achieve the windmill regime by tuning the current exactly 
to the "marginal" value such that C(/^) = 0. Note 
that this value corresponds to a spin transfer strength 
of Uw ~ a Jm 17712, and since a <^ 1 the required spin 
torque is much smaller than the exchange torque. The 
situation is similar to the switching regime, where spin 
transfer effect works against the magnetic anisotropy. In 
both cases critical values of spin torque are proportional 
to the small Gilbert damping coefficient. 

The original discussion of the windmill regime in Ref. [2 
assumed J = 0, ai_2 = 0, and mi = m2- It was found 
that magnetic moments rotate in the plane spanned by 
vectors mi and m2 at the initial moment, and the angle 
9 between them remains constant. How will the windmill 
motion look in the general situation? At C = the total 
magnetic moment is conserved, = m^ -|- ml + 2 (mi • 
m2) = const, thus 6 is constant in general case as well. 
Since magnitudes of mi and m2 are also fixed, constant 
6 implies that both vectors will rotate with the same 
angular velocity. 



(11) 



(cf. the theorem on the motion of a rigid body with a 
fixed point) . To find a; we expand it in the basis of vectors 
(i>',mi,m2) as a; = ai^* + feimi + 621112 with unknown 
coefficients a and 61^2- Substituting this form of a; into 
Eqs. (dU, using expressions ([7]) for rii^, the fact that for 
the marginal value of current one has Bi = —B2 ^ B^, 
and properties [u x mi] ~ — li, [1/ x m2] = I2, we find 
a = —Bu], bi = A2, and 62 — Ai 



u) = Byj[mi X m2] -t- A2I111 + ^im2 



(12) 



Since Cj ~ [a; x a;] = 0, u; is an invariant of motion, 
determined by the initial conditions. 

Since ai_2 1 and Uw olJ at the marginal point, 
we can make approximations in expressions ([S]) and use 
Ai w A2 K, 7J, Bi K, j{aiJmim2 + Ui)/m\m2, B2 ~ 
j{a2Jmim2 — U2)/mlmi. Equation C = Bi + B2 — 
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then gives the following values at the marginal point 



with 
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mia2 + 'm2ai 
mi52 - m2gi 



Jmim2 



(13) 
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ai92 + (^291 

7711(72 - "7,251 
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Note that approximation Uw ^ aJ violated when pa- 
rameters are close to the degenerate situation 7711^2 — 
= 0. This is the situation when C is independent of 
the current and the windmill regime cannot be achieved. 
Far away from the degenerate situation one has 



7 J 



aiff2 + Qi2gi 

7)7.132 — "^251 ' 



[mi X m2] + (mi + m2) 



(14) 

The first term in parentheses is smaller than the sec- 
ond one by a factor of a ^ 1. Thus vectors m.^ pre- 
cess approximately around the total magnetic moment 
M = mi -|- m2. In the presence of damping and spin- 
transfer M is not conserved and performs a small angle 
rotation around the constant vector lj. In this respect 
the motion is very different from Slonczewski's situation 
at J = 0, where M was performing 360° rotations around 
[mi X m2]. The J = rotations, however, require the ful- 
fillment of a condition gi/[(l -I- 0:^)7711] = g2/[{l + a\)m2], 
which is, in particular, satisfied in a completely symmet- 
ric valve considered in Ref. 0. 

Finally, we return to Eq. ^ and investigate the C 7^ 
case. It is convenient to rewrite ^ in terms of a; = cos 6 



X = Cmim2{l — x'^) 



The solution reads 



x{t) = cosd{t) — tanh 



sgn[C] 



\C{I)\mim2 



(15) 



and parameter tg determined by the initial angle, cos — 
tanh(to/r*). We conclude that as t ^ 00 the system ap- 
proaches a coUinear configuration with a current depen- 
dent characteristic time r*(J). The latter diverges in the 
vicinity of the marginal current Iw 

In conclusion, we studied the motion of a two layer 
spin-transfer device with zero magnetic anisotropy. We 
show that in the presence of damping, layer asymmetry, 
and exchange interaction between the layers the windmill 
rotation decays with characteristic time constant . The 
decay time depends on the current pumped through the 
device and diverges at a "marginal" current Iw For / ^ 
Iw the system reaches either a parallel or an antiparallel 
state after a transient period. Exactly at the marginal 
point / = In, the system performs a perpetual generalized 
windmill motion. 

Interestingly, precession motion analogous to the wind- 
mill regime was also found in multilayers and bilayers 
with magnetic anisotropyii^ In those systems it exists 
not at a singular point, but in the whole range of cur- 
rent values. Thus, rather unexpectedly, anisotropy can 
be advantageous for the windmill regime. 

Finally, coming to the discussion of the current in- 
duced ferromagnctism, we see that in a two magnet de- 
vice current can induce both ferromagnetic and antifer- 
romagnetic order. However, the situation with only two 
magnets can be special, and it is necessary to consider 
devices with three and more magnets to predict what 
happens in the system of many isotropic paramagnetic 
impurities under the influence of spin-transfer torques. 
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